On optical black holes in moving dielectrics 
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We study the optical paths of the hght rays propagating inside a nonlinear moving dielectric 
media. For the rapidly moving dielectrics we show the existence of a distinguished surface which 
resembles, as far as the light propagation is concerned, the event horizon of a black hole. Our 
analysis clarifies the physical conditions under which electromagnetic analogues of the gravitational 
black holes can eventually be obtained in laboratory. 
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One of the most striking and well established results of 
general relativity theory is the prediction of a bounded 
region in space from where nothing can escape, not even 
light. Such an object, a black hole, arises as the spher- 
ically symmetric solution of Einstein equations, and its 
analysis on the basis of the semi-classical quantum field 
theory was considered to be fundamental for the under- 
standing of quantum gravity. In particular, in the con- 
text of the semi-classical gravity. Hawking suggested that 
black holes can evaporate Since the temperature of 
the Hawking radiation is very low, astrophysical events 
are unlikely to give an evidence for such a phenomenon. 

During the last decades efforts have been put into a 
development of the alternative physical models in which 
such a prediction could be experimentally verified under 
laboratory conditions. This class of models was called 
'analog gravity' 0, and it deals with _the_ propagation of 



and with the 
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sound waves in moving fluids ||3|, U, 
electrodynamics of nonlinear media 
among others (see for instance Ref. 13| and references 
therein). 

The search for electromagnetic black holes represents 
an interesting topic in the broad area of analogue grav- 
ity, but the effective construction of such structures still 
remains an open problem. In a recent work Jl^ a certain 
progress has been achieved in this direction. In particu- 
lar, it was shown that the birefringence effect occurring 
in moving nonlinear dielectrics can be described in terms 
of the optical geometries (also known as the electromag- 
netic effective geometries). Accordingly, the propagation 
of the electromagnetic light rays is governed by the opti- 
cal metric tensors: 

gii" ^ + (e^ - l)u^'u'' - fiu^^'e''^°'F^pu'^ 

_i^^e"('^^'')%, (1) 
^ rJ^^'^ + (e//- l)u^u^ (2) 

Here 77'"' is the background (flat) spacetime metric, and 



the magnetic permeability fi is assumed to be a constant. 
The derivative e''" = de/dFfj^^ of the dielectric permittiv- 
ity e with respect to the electromagnetic field measures 
the nonlinear effects. The motion of the medium is de- 
scribed by the velocity vector measured in units of the 
velocity of light in vacuum (i.e., c = 1 except specified 
otherwise). The meaning of the Eqs. (0),(§) is that the 
wave vectors K\ — d\Ti lie on the null geodesies of the 
above two optical metrics, depending on their polariza- 
tion state. Here S denotes the corresponding wave-front 
surface. 

The aim of the present work is to show explicitly how 
the analogue event horizons can be built from the above 
metric structures. We analyze the propagation of the two 
light rays described by the optical metrics (|l|),(|^). De- 
spite of birefringence, the horizons for both polarization 
modes are shown to coincide. The spherical configura- 
tion is discussed, and we find an analogue event horizon 
resembling the Schwarzschild event horizon. Further, we 
study the cylindrically symmetric configuration as a pos- 
sibly simpler structure to be built in a laboratory. 

Let us now examine the light propagation inside a 
moving nonlinear dielectric matter. Consider at first the 
spherically symmetric case with the background metric 
ds^ = dt^ — dr'^ — r'^dO'^ — sin^ 9 dip^ in which the di- 
electric flows radially with the velocity = 7 (1, 0, 0), 
where 7 = (1 — /3^)~^/^. We assume a point charge fixed 
at the origin, and correspondingly the only non-zero com- 
ponent of the electromagnetic field is then = E,.. 
Thus e°iFoi = E{de/dE) = e', with E = \Er\. The 
optical metric ([T]) then reads 
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Note that g^" — g'^'^\e'=o- For spherically symmetric 
flows /3 — l3{r), the two metrics g'^ evidently admit the 
two Killing vector fields ^ = {d/dt,d/d(p}. Written in 
components, — Sj^ with a ~ 0,3. The light then prop- 
agates along the null geodesies with the tangent vectors 
j^ii = — {t,f,9,if) in such a way that £,qK^ — Kq 
and ^ijKfi = K3 are constants along the corresponding 
optical path. Here the 'dot' denotes the derivative with 
respect to the aSine parameters s along the optically null 
geodesic lines. For simplicity, we choose the radial prop- 
agation with K2 = = K3. The path of the light ray is 
then determined hy t — r{nzL f3) / {nP ±1) . Here n = ^/Jle 
is the refraction index of the medium in the absence of 
external fields. The plus (minus) sign corresponds to the 
wave propagating from (to) the origin. The radial coor- 
dinate of the photon then changes in time as 



dr n(3 ± 1 
'di ~ n±(3 ' 



(7) 



For the case of the matter flow directed outwards, i.e., 
/3 > 0, then the incoming light (with r < 0) apparently 
"freezes" at the radius rh such that nf]\r^ = 1. The same 
happens for an outgoing light in the inwardly directed 
flow: the "freezing" occurs then at n(3\r^ = —1. Both 
cases are compactly summarized in 



1. 
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It is worthwhile to note that "freezing" is absent for the 
light propagating in the same direction as the moving 
dielectric medium. 

There is a clear physical interpretation of the above 
result. Recall that P — v/c tor the 3- velocity v of the 
matter flow, and Vc = c/n is the velocity of light in a 
medium with refraction index n in the absence of external 
flelds. Then we can straightforwardly recast the equation 
(0) into the form of the relativistic transformation of the 
velocity: 



dr V ±Vc 
dt l±vvc/c^ 



(9) 



The 'freezing' of light takes place at a surface on which 
the velocity of matter v becomes equal to the velocity 
of light Vc in the medium. The geometrical interpre- 
tation is also clear: we find that the metric coefficient 
gl}{rii) = vanishes at the radius r^, see Eq. (|^). The 
same conclusion is true also for the second optical metric, 
gi^{rh) = 0, i.e., both polarization modes (|l|) and (||) be- 
have similarly. For the radial motion considered above, 
the stationary metrics given in the Eqs. (|l|), (||) can be 
rewritten into the static "black hole" form 



dsl = , 
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with = 7^(1 — n^/3^)/^(e + e'), by means of the co- 
ordinate transformations 



9± 

dt = dt -\ dr, 

9± 



df = drJ{gfy-g^^g}^ 



(11) 
(12) 



We note, for the case here described, that 5^ is propor- 
tional to g}^ . 

From the above analysis, which generalizes the pre- 
vious calculations ijlTI , we find that either an ultra- 
relativistic motion of the matter or a highly refringent 
medium [|l5|, |l^, as it is found in Bose-Einstein con- 
densates (see, however, a relevant discussion in Ref. |^), 
is needed in order to display such a horizon structure. 
Highly refringent media are usually quite dispersive, and 
thus the effective horizon would occur in this case only 
for a narrow range of frequencies. A similar feature oc- 
curs also for sonic analog models, where the horizon is 
defined only in a narrow range of low frequencies. A 
practical realization of the analogue black hole configu- 
ration requires that either n = n{r) or /3 = /3(r), or even 
both. 

Recently, it has been proposed that Hawking tem- 
perature is a purely kinematic effect that is generic to 
Lorentzian geometries containing event horizons , and 
thus being dependent only on the effective metric struc- 
ture. For the spherically symmetric solution the Hawking 
temperature T of the analogue black hole is straightfor- 
wardly calculated for the metric (§)-(§) as 



T : 
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Equation ( p^ ) shows that the nonlinearity may signifi- 
cantly contribute to this temperature. It should be noted, 
however, that the complete understanding of the physical 
meaning of the Hawking temperature associated with an 
analogue black hole is still an open question Q , and it can 
be settled only after a detailed analysis of the radiation 
processes at the horizon. It should be stressed that the 
approximation of geometrical optics becomes unreliable 
for the discussion of modes whose wavelengths are com- 
parable with the size of the horizon, as occurs in Hawking 
radiation processes. 

Experimentally though, it appears to be a rather diffi- 
cult task to maintain a stationary spherically symmetric 
and inhomogeneous flow. In order to exhibit a more re- 
alistic configuration, we will now focus on the cylindrical 
symmetry. A particular case of such a configuration is the 
vortex matter fiow which was discussed recently [^8[ ^ . 
Let us consider a more general situation of a rotating 
dielectric body subjected to an electric field directed 
along the axis z of rotation |2^. The background met- 
ric is then of the form ds^ — dt^ — dp^ — p^dLp^ — dz^ . 
The electromagnetic field has the only non-zero compo- 
nent = and wc denote e^^Fos = E{de/dE) = e'. 
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where E — \Ez\. For the matter 4- velocity we have 
= 7(l,/3,w,0), with 7 = 1/^1-/32 - p^cj^, while lu 
and P being arbitrary functions of the radial coordinate 
r. The optical geometry (|^) now reads 
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For the motion under consideration, the above metric 
components, Eqs. (14)-(|Tg|), depend only on the radial 
coordinate p. Then, for this geometry we have the 
three Killing vectors d / dt^d / dip,d / dz. In components, 

= '^ith a = 0,2,3. The null geodesic equations 
K^^\/^Kx = for the wave vector K\ in the optical met- 
ric g'^'^ are then more easily solved by using the fact that 
^'^Kfi — const for each Killing vector That is, the 
covariant components Kq, K2, K3 of the wave vectors 

— i^^ — {t, p, ip, z) are constant along the correspond- 
ing geodesic lines. The remaining component Ki is ob- 
tained from the null condition K^j^K^ = 0, and the 
geodesic equation for the radial coordinate reads 



^33 



e 



(19) 



p = ±J{gl^Ko + gl'K^r - 9}}[2g^3K^K2 + g'_'{K,y + gl^K^y + g^.^K^n 



(20) 



where the plus (minus) sign corresponds to a light ray propagating from (to) the symmetry axis z. Also, z = g_ -fCa 
and 



t = 



.9" 



(21) 
(22) 



By specifying the dynamics of the matter, i.e. the func- 
tions (3 and w, all the metric components g^ will be 
given functions of p. Provided p = p{s) is obtained from 
Eq. (|2^), the three other coordinates are then found by 
quadrature. 

Let us now analyze the ordinary ray described by g'^'^ , 
which is obtained from the above computations when 
we set e' = 0. For simplicity, we also consider the pla- 
nar initial conditions for the light propagation such that 
K2=0 = K3. Equations (|20|) and (^ then yield 



dp 
'dt 



(n2 - 1)7/3 =F ^n2(l - /32) _ ^2^2 



(23) 



As compared to the spherical case, Eq. (|7|), the solution 



of Eq. (23) displays the same qualitative behavior for 
the cylindrical configuration. The cylindrical horizon is 
located at p = ph defined from 



p2^2 



0. 



(24) 



Purely vortical motion of the fluid (with P — 0) does not 
produce an analogue event horizon ]l^ , since ph from 
Eq. (|2^) would lie beyond the limit of applicability of a 
rigid body in special relativity. 



We note that dip /dp diverges at the surface defined by 
Eq. (p^). As a result, the incoming/outgoing geodesic 
light rays spiral towards a horizon radius ph, which can 
be explicitly demonstrated by the numerical integration 
of dip /dp, see Fig. |^. Moreover, the condition of a hori- 
zon for the optical metric g'/i'^ is precisely given by the 
Eq. (|^. In terms of the optical metrics, Eqs. (|l^-(^, 
the position of the horizon is defined by the equation 
gj}{ph) = 0. In Fig. |l|, we also depict the extraordinary 
ray which presents a qualitative behavior quite similar to 
its ordinary counterpart. Although the calculations look 
somewhat more involved when the above initial condi- 
tions are not fulfilled, it can be shown that the horizon 
condition ( p4[ ) remains unchanged. Therefore, we con- 
clude that, for a given stationary flow configuration, the 
analogue horizon structure is of geometrical nature, as 
soon as it depends neither on the initial conditions nor 
on the polarization of the propagating light rays. 

Summarizing, Eq. (|2j) demonstrates how an experi- 
mental realization of a dielectric analogue horizon might 
be understood in terms of the parameters n, /?, lu which 
describe the dielectric and kinematic properties of the 
medium. We expect that such horizons can be observed 
for stationary inhomogeneous kinematic configurations of 



4 




FIG. 1: The light path for both polarization modes are shown for the cyhndrical ('black string') configuration. The left/right 
figure corresponds to a dielectric fiuid fiowing outwards/inwards, with radial velocity /3out(p) = / p and /3in(p) — — lO^^p, 
respectively. The solid lines depict the ordinary rays, whereas the dashed lines represents the extraordinary rays. For simplicity, 
we assume in the plots a rigid rotation with ui — 10~^, and the refraction index is n = 10^. Besides, we take ^te' = 10^" for the 
extraordinary rays. (The choice of a large value for e is but illustrative, since the typical values would produce the rays too 
close to be distinguished from one another in the plots.) The thicker lines correspond to the analogue horizons in both cases. 



the dielectric matter flow. The cylindrical configuration 
is favored due to the combined effect of both the refrac- 
tion index n and also of the vorticity lo of the medium 
which yields a smaller threshold for the radial velocity 
j3. The horizon structures in both cases are shown to 
be independent of the presence of nonlinearities in the 
permittivity tensor. 

As a final remark, it is worthwhile to mention that the 
present scheme may provide room also for the formation 
of sonic horizons (dumb holes [|| ) for the cases in which 
ultra-sonic velocities are achieved. 

In the present paper, the study of light propagation in 
a nonlinear moving dielectric medium has revealed the 
existence of a surface which displays some of the proper- 
ties of the event horizon of a black hole. We thus confirm 
the possibility of creating electromagnetic analogues of 
gravitational black holes in laboratory conditions. 
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